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Abstract 

We study the stability of solution branches for the Lichnerowicz-York equation at mo- 
ment of time symmetry with constant unsealed energy density. We prove that the weak- 
field lower branch of solutions is stable whilst the upper branch of strong field solutions 
is unstable. The existence of unstable solutions is interesting since a theorem by Sattinger 
proves that the sub-super solution monotone convergence method only gives stable solu- 
tions. 



1 Introduction 

The loss of uniqueness is a salient feature of many nonlinear systems and is often accompa- 
nied by a loss or "exchange" of stability as two solution branches meet. A number of recent 
papers have looked at the uniqueness of solutions of conformal formulations of the Einstein 
constraints. 

Recent numerical work by Hoist and Kungurtsev in BH has confirmed the assumption made in 
lfT9ll that the linearisation of the (constant) energy density moment of time symmetry Lichnerowicz- 
York equation has a one dimensional kernel and that a quadratic fold results giving two solu- 
tions for subcritical energy density p < p^, one solution at the critical value p = Pc and none 
for energy density greater than a critical value p > Pc- 

In [fT3l. Maxwell studied the conformal formulations of the Einstein constraint equations with- 
out the usual constant mean curvature condition (CMC) so that the constraints do not decouple. 
One motivation for that study was to examine the uniqueness of solutions in the far-from- 
constant mean curvature regime now that existence results for this case exist (see [H, [[T4l ). 
He found interesting non-existence and non-uniqueness results showing that the constraints 
are ill-posed beyond the non-uniqueness that is introduced when one couples the lapse fixing 
equation to the four constraint equations, as in the extended conformal thin sandwich (XCTS) 
formulation (see ifTSl . flU, [[T9l ) and some constrained evolution schemes (see [|T6l . [|6l for a 
resolution of this scaling problem). 

The inherent ill-posedness of conformal formulations of the constraints found in [fT3l is worthy 
of further analysis. In this work we continue our bifurcation analysis of the Einstein constraints, 
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begun in [[T9l with an analysis of the XCTS system, by studying a familiar problem and intro- 
ducing the important related issue of the stability of the solutions obtained. The ultimate goal 
of this program is to better understand the mapping between general free initial data and solu- 
tions of the constraint equations in their various conformal formulations. See also the recent 
work of Hoist and Meier [ fTOl in this regard. 

In section 2 we study the solution branches of the constant density star found in f2\, f[9\ and 
m. We prove that the lower branch of weak- field solutions is stable whereas the upper branch 
is unstable. (W also prove that the kernel of the linearisation is one-dimensional, which was 
shown numerically in f9l). In section 3 we follow [fT2l in applying Liapunov-Schmidt meth- 
ods, as in [T9], to the stability analysis and find that the exchange of stability that occurs at the 
critical point of a fold is in fact generic under mild non-degeneracy conditions. The existence 
of unstable solutions to the constraint equations is interesting in itself because the most popu- 
lar method for proving existence, the sub-super solution monotone convergence method, only 
yields stable solutions (a result proven by Sattinger in ifTSll ). 



2 Stability of solutions 

We shall be concerned with the Hamiltonian constraint, the Lichnerowicz-York equation, at 
moment of time symmetry, for a conformally-flat background metric and an unsealed constant 
energy density p: 

F(0,p) := V20 + 27rp0^ = (2.1) 

with > 0. We work in spherical symmetry in this section so that = + f ^- For 
simplicity we assume that p is zero outside a ball of radius r = 1 with boundary conditions 
0(1) = l,andf (0) = Oas in L9J. 

We are interested in the stability of stationary solutions of the following parabolic problem: 

— = VV + 27rpz/^ (2.2) 

with iy{t = 0, x) = (f)o{x). 

We are interested in whether, for initial data (f)o{x) close to the stationary solution, the solution 
to this parabolic problem tends to the stationary solution as t — > oo. We take 

p = + e9i 

where z> is the stationary solution and 9i is the eigenfunction corresponding to the principal 
eigenvalue pi (the smallest eigenvalue) and e is a small constant. Then substituting this into 



(2.2) gives 

— Bf) 

^ V^Oi + 107rp(j)% = -piOi 



dt 

to first order in e. 



Definition We will say that a stationary solution to ( |2.2| ) is stable if the principal eigenvalue of 
the linearisation 



+ 107rp(j)% = -pkOk, (2.3) 
with 9k{l) = and 6*^(0) = satisfies /ii(V^ + lOvrp^^) > and otherwise that it is unstable. 
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In the next Lemma we collect some facts about the eigenfunctions and eigenvalues of ( |2.3| ). 
Lemma 2.1 



1. All the eigenfunctions of (|2.3|) are orthogonal and all eigenfunctions except 9i have 



nodes. Furthermore, the principal eigenvalue is simple i.e. iff is any solution of 

VV + lOvrp^V = -/ii/, (2.4) 

then f is proportional to 9i 

2. The principal eigenvalue has a variational characterisation given by the Rayleigh quo- 
tient: 

/ii(V + a{x)) = mmrj^o^D J^^fd^ ^ -* 

where D is the set of smooth functions satisfying the boundary conditions r]{l) = 0, 

S(o) = o 

3. Ifa{x) > b{x) with a{x) > h{x) in a subset of positive measure then 

/ii(V^ + a{x)) < /ii(V^ + h{x)) 



4. Similarly, ifa{x) < b{x) with a{x) < b{x) in a subset of positive measure then 

/xi(V2 + a{x)) > pi(V2 + b{x)). 



Proof 1. It is easy to show that the eigenfunctions of (2.3 1 are orthogonal: 



/ 9k9jdV= [ {9kiV% + 107rp(t>%)~9jiW% + 10TTp(j)%))dV 
Jv Jv 



{9kV9j ~9jV9k).ndS 

dV 

= 

and we have used Green's second identity with n the outward pointing normal to the surface 
element dS at r = 1. We know that the eigenfunction corresponding to the first eigenvalue fii 
is of definite sign so may be taken to be positive (see[7J). The orthogonality relation above then 
says that all higher eigenfunctions must have nodes. The proof that /ii is simple may be found 
in [|71 for example. 

2. This well known result may be found for example in ([0). 

3. and 4. follow directly from the above result. 
QED 

In the authors worked in an unbounded domain which complicates the analysis of the spec- 
trum of the linearisation. We work in a bounded domain corresponding to a ball of radius r = 1 
with 0(1) = 1. It is straightforward, using the ideas in [|2l, to generate a one parameter family 
of solutions with this boundary condition: 
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where the boundary condition gives 

3 / a \^ 
ZTT y 1 + « / 

Note that the maximum value of p := occurs at = 1 and that the limit a — )■ oo corresponds 
to a lower branch of solutions and the limit a — )• corresponds to an upper branch of solutions. 

It is easy to check that 

(^k{r) = (2.8) 
(1 + r^)2 



satisfies (2.3) with a = 1 in p and = 0. This eigenfunction has no nodes and therefore 
corresponds to the principle eigenvalue found numerically in [9|. So we have that /ii = 0, we 
have found the principal eigenfunction at the critical value pc, = \. Since this eigenvalue 
is simple we have that the kernel of the linearisation at the critical energy density pc is one 
dimensional. 

In [|2l, and above, global branches of solutions were found corresponding to a quadratic 
fold. For each p < pc we have an upper solution 0{7 and a lower solution 0^. (As in [2J, a 
simple calculation reveals that there are apparent horizon's on the upper branch of solutions at 
the location r = a, and none on the lower branch). 

Before showing that the upper solutions (\)u are unstable and the lower solutions 0^ are stable, 
we motivate our analysis by the following observation. 



If we multiply the linearisation (2.3) by m := — 1 and integrate twice by parts we get the 



following identity, having used (2.1 ): 

27r y (p0^^i(4m - 1)) dv = -f^ij ^ 



ludv. 



With 9i the first eigenfunction and therefore positive and m > since > 1 we see that for 
weak-field solutions (4m — 1) is likely less than zero so that /ii > and the solution is stable. 
But we can imagine a strong field solution whereby (4m — 1) > so that < and we have 
an unstable solution. 

Theorem 2.2 

1. The lower branch of solutions, i.e. (pir; a) with 1 < a < oo, is stable. 

2. The upper branch of solutions, i.e. 0(r; a) with Q < a < 1, is unstable. 

Proof 1. We know that pi(V^ + lOvrpc^^) = 0. Motivated by this we examine the difference 



in the "potentials" in the linearisation (2.3) 



p(a)0(a) - p.0, = - - J = ^ (^(,2 + i)2(,2 + «2)2 J (2.9) 
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Note that r G [0, 1]. On the lower branch of solutions 1 < a < oo, so that p(a)0(a)'^ — pc0c < 
0, so Lemma 2.1.3 implies that 



/ii(V^ + lO7rp(a)0(a)^) > /^(V^ + 107rp,(/)^) = (2.10) 
so the lower branch of solutions is stable. 

2. We now proceed to show that the upper branch of solutions is unstable. The comparison test 
above yields an indefinite result so we cannot implement Lemma 2.1.4. 

Instead we utilise the fact that 

Pi(V + lOn p(f) ) = Tmrirj^oeD- J^rfdv 

= mzn^^oeD J rj'^dv ' (2-12) 

Using a trial function rjT E D then provides an upper bound for the principal eigenvalue corre- 
sponding to a particular choice of < a < 1, i.e. pi(a) < R{a) where R{a) is the Rayleigh 
quotient above. 

We shall choose 

l-r2 

riT = w (2.13) 

and use MATLAB to perform the integration exactly. This yields a rather complicated expres- 
sion given in Appendix A. A graph of R{a) is shown in Figure [T]and some explicit values for 
R{a) are given in Table 1 in Appendix A. We clearly see that for < a < 1, the upper branch 
of solutions, that fii{aupper) < R{ceupper) < so we have that the upper branch of solutions is 
unstable. 

QED 



3 Stability of solutions via Liapunov-Schmidt reduction 

Liapunov-Schmidt (LS) methods were applied to the Lichnerowicz-York equation with an un- 
sealed source in [19] and more recently in [fTOl. Here we follow [12] and show that locally the 
exchange of stability witnessed in the previous section is in fact generic for fold-type branches. 

We again consider a nonlinear equation F(x, p) = with x E X, a Banach space and p G M, 
(F : X * M — 7- Z) with X C Z and assume that zero is a simple eigenvalue of DxF{xo, po). 
We let 9 denote the one-dimensional kernel and denote by 6* the cokemel of D^F^Xq, po) 
such that J 66* = 1 (we assume the linearisation has a Fredholm index of zero). The LS 
construction allows us to construct a continuously differentiable curve of solutions (x(s), p(s)), 
with s G {-S, S) and (a;(0), p(0)) = [xq, po), such that F{x{s),p{s)) = for all s G {-6, 6). 
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The Rayleigh quotient upper bound on n for fixed a 




a 



Figure 1 : The Rayleigh quotient as a function of a for the trial function r]T 

For the LS projection operators and splittings of the range and domain we follow the notation 
of [fT2ll section 1.4-1.7, which is equivalent to [fT9l . In particular we have that 

Z = RiD^Fixo, po)) © NiD^Fixo, po)) (3.14) 
X = N ®{RnX). (3.15) 

Note that if the following nondegeneracy conditions hold: 

D,,F{xo,Po)[0,9] i RiD,Fixo,po)) (3.16) 
DpF{xo, Po) ^ RiD^F{xo, Po)) (3.17) 

then it is proven in [[T2l that 

p(0)=0, p(0)^0 (3.18) 

where a dot denotes differentiation with respect to s, and so the tangent vector to the solution 
curve (s(s), p(s)) at (xq, po) is (9, 0) (see corollary 1.4.2 in [[El) and we have a turning point 
or fold (sometimes called a 'saddle-node'). 

To analyse the stability of the 'branches' of the fold (upper and lower) we note the following 
result (Proposition 1.7.2 in [fT2l ): 

Theorem 3.1 

There exists a continuously dijferentiable curve of perturbed eigenvalues {p(s) |s G (—5, 5) p(0) = 
0} m M such that 

D,F{xis), p{s))ie + wis)) = pis)ie + wis)), 

where {w(s)|s G i—6,6), w(0) = 0} C R H X is continuously dijferentiable (and the size 
of the interval (—5, 6) is not necessarily the same as in the solution curve above, but possible 
shrunk). In this sense, p(s) is the perturbation of the critical zero eigenvalue ofD^Fixo, po). 
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(Note that we used the opposite sign convection for the eigenvalues of the linearisation in 
section 3 where D^F^xq, pq)9 = —fi9). 



It is then straightforward to show that (see 1.7.30 in [fT2ll \ 

m = J D^,F{xo,po)[0,e]e* ^ (3.19) 



and 



m = - (/ DpF{xo,po)e*^ p(0). (3.20) 



Now /i(0) = 0, /i(0) 7^ implies that ji(s) changes sign at s = so that the stability of the 
solution curve changes at the turning point. 

Applying this to the constant density problem of section 3 we see that for small s E (—5, 5) 
(and with a and 7 two known constants obtained from ( |3.19[ ) and ( |3.20[ )) 



p(s) =Pe-7V + 0(s3) 

which agrees with [[T9ll and section 3 above, and 

H(s) = a^s + O(s^) 

which tells us that the lower branch of solutions s G (—5, 0) is stable and the upper branch 
s G (0, 5) is unstable, as expected. 

We conclude this section by mentioning a limitation on the application of the sub-super so- 
lution method that this work reveals. Solutions to elliptic equations obtained via the sub- 
super solution method are stable in the sense described in section 2 (see [Ol, ifTSl ) so that 
the upper branch of solutions studied here are unattainable by this method. The general CMC 
Lichnerowicz- York equation with an unsealed source on an asymptotically Euclidean manifold 

- r(/) + a0- V 27rp05 = (3.21) 

was studied in flU using the sub-super solution method, where r is proportional to the scalar 
curvature and a the traceless part of the conformally transformed extrinsic curvature squared 
(see section XII of that work). They found an open set of values of a and p such that existence 
could be proven-uniqueness was not proven. 

If, as seems likely from this work and [[TSl . O, [[T9l .[l9ll and [[TOll . that upper and lower branches 
of solutions exist for this equation for some combination of a and p then the sub- super solution 
and monotone convergence method will only converge to a stable solution. So if an upper 
branch of solutions is unstable, as above, then the sub-super solution method will not yield 
it. These comments are also relevant to the Einstein-Scalar field CMC Lichnerowicz- York 
equation as studied in [[5]| (see Theorem 8.8 in that work). 



4 Discussion 

In this work we have concentrated on the non-standard case of an unsealed fluid with no mo- 
mentum. (For a complete discussion of the role of conformal scaling of fluid sources see [|3]| or 
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IH). It is not of purely academic interest however because it serves as an excellent model of 
a poorly scaled system such as the XCTS system and also warns of the dangers of not scaling 
the extrinsic curvature as was common in numerical evolutions which started from moment of 
time symmetry initial data, see section 5 in [|19J . 



In studying (2.1 1, we have been looking at the geometric problem of finding a conformal factor 
that maps from a scalar flat metric to one with scalar curvature equal to 16np (for more details 
see the classic paper [[TT| where the authors allow a non-zero scalar curvature background 
metric). This is closely related to the Yamabe problem of finding a conformal factor that maps a 
given metric to one with constant scalar curvature. This has yielded many insights, particularly 
in relation to the existence theory of solutions to the Lichnerowicz-York equation, see [Qj for 
an excellent review. 

It is clearly important to understand the strengths and limitations of solution methods. With 
this in mind, future work will determine the stability of non-unique solutions found in other 
works, for example in [fTSl . [fT3l. 

The Lichnerowicz-York equation with an unsealed source also has a variational formulation 



since ( |2.1| ) is the Euler-Lagrange equation obtained by varying the following functional: 



dv (4.22) 



It would be interesting to check which solutions (stable/lower or unstable/upper or both) apph- 
cation of the Mountain Pass theorem to this toy model yields. 
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Appendix A 

We used MATLAB to perform the following integration for the Rayleigh quotient for the trial 



function riT given by (2.13 1 



fii{a)<R{a) (4.23) 



(4.25) 



{3tan^^{l/a)a^ — So;'^ + — 6tan^^(l/a)a^ — + 3a + 3tan~^(l/a)) 
((a^ + 2a2 + I)a2(15tara-i(l/a)a4 - ISa^ + Qtan-l{l/ a)a^ -a- tan-\l/a))) 

(4.26) 



8 



Some specific values of the Rayleigh quotient are given in Table 1 below. 



Table 1: 



a 


Rayleigh quotient 


10,000 


10.500 


5 


9.804 


1 





0.99 


-0.138 


0.75 


-4.548 


0.5 


-14.057 


0.25 


-53.713 


0.01 


-3.00x10^ 



Using MATLAB again we can confirm that the Rayleigh quotient horizontally asymptotes to a 
value of 10.5 as a — > oo. 
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